Measurement of the quantum zero-point motion of a nanomechanical resonator 
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We present optical sideband spectroscopy measurements of a mesoscopic mechanical oscillator cooled near 
its quantum ground state. The mechanical oscillator, corresponding to a 3.99 GHz acoustic mode of a patterned 
silicon nanobeam, is coupled via radiation pressure to a pair of co-localized 200 THz optical modes. The 
mechanical mode is cooled close to its quantum ground state from a bath temperature of T/, 18 K using 
radiation pressure back-action stemming from the optical pumping of one of the optical cavity resonances. An 
optical probe beam, resonant with the second optical cavity resonance, is used to transduce the mechanical 
motion and determine the phonon occupancy of the mechanical mode. Measurement of the asymmetry between 
up-converted and down-converted photons of the probe beam yields directly the displacement noise power 
associated with the quantum zero-point motion of the mechanical oscillator, and provides an absolute calibration 
of the average phonon occupancy of the mechanical mode. 



The Heisenberg uncertainty principle, one of the fundamen- 
tal consequences of quantum theory, restricts the certainty with 
which the position and momentum of an object may be simul- 
taneously known and defined. Consequently, an object con- 
fined to a local potential possesses a non-zero ground state en- 
ergy associated with random quantum fluctuations of its posi- 
tion. For mechanical systems of our daily experience, this so- 
called quantum zero-point motion is masked by thermal noise 
resulting from interaction with the environment. To observe 
quantum motion, the dual and antagonistic feats of isolating a 
mechanical system from its environment whilst measuring its 
position must be achieved 1 1 -3|. Recently, coupled electro- 
and opto-mechanical systems have been cooled into their me- 
chanical quantum ground state by a combination of cryogenic 
pre-cooling and radiation pressure back-action ||4l|5l. In this 
work, we extend these results by laser cooling the acoustic 
mode of a silicon optomechanical crystal resonator to near its 
quantum ground state (occupancy 2.6 ± 0.2 phonons) while 
simultaneously monitoring its motion with a probe laser res- 
onant with a secondary optical mode of the cavity. Operat- 
ing in the resolved sideband regime, the spectral selectivity 
of the secondary cavity mode is used to separate, measure, and 
compare the intensity of the blue- (anti-Stokes) and red-shifted 
(Stokes) sidebands created by the mechanical resonator's mo- 
tion. The observed asymmetry in the probe sideband am- 
plitudes provides a direct measure of the displacement noise 
power associated with quantum zero-point fluctuations and al- 
lows for an intrinsic calibration of the phonon occupation num- 
ber of the nanomechanical resonator. 

Experiments with trapped atomic ions and neutral atoms (6]- 
[8l, dating back several decades, utilized techniques such as re- 
solved sideband laser cooling and motional sideband absorp- 
tion and fluorescence spectroscopy to cool and measure a sin- 
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gle trapped particle in its vibrational quantum ground state. 
These experiments generated significant interest in the coher- 
ent control of motion and the quantum optics of trapped atoms 
and ions|9|, and were important stepping stones towards the 
development of ion-trap based quantum computing ITOl [TTIl . 
Larger scale mechanical objects, such as fabricated nanome- 
chanical resonators, have only recently been cooled close to 
their quantum mechanical ground state of motion (H O \V2\ - 
17 1 . In a pioneering experiment by O'Connell, et al. |[T6ll , 
a piezoelectric nanomechanical resonator has been cryogeni- 
cally cooled (7], ~ 25 mK) to its vibrational ground state and 
strongly coupled to a superconducting circuit qubit allowing 
for quantum state preparation and read-out of the mechan- 
ics. An alternate line of research has been pursued in cir- 
cuit and cavity optomechanics where the position of a 
mechanical oscillator is coupled to the frequency of a high- 
Q electromagnetic resonance allowing for back-action cool- 
ing [18, J9| and continuous position read-out of the oscil- 
lator. Such optomechanical resonators have long been pur- 
sued as quantum-limited sensors of weak classical forces C]- 
[3l[T4l[20J, with more recent studies exploring optomechanical 
systems as quantum optical memories and amplifiers EDEH, 
quantum nonlinear dynamical elements | 25 1, and quantum in- 
terfaces in hybrid quantum systems |26-29|. 

Despite the major advances in circuit and cavity optome- 
chanical systems made in the last few years, all experiments 
to date involving the cooling of mesoscopic mechanical os- 
cillators have relied on careful measurement and calibration 
of the motionally scattered light to obtain the average phonon 
occupancy of the oscillator, {n) . Approach towards the quan- 
tum ground state in such experiments is manifest only as a 
weaker measured signal, with no evident demarcation between 
the classical and quantum regimes of the oscillator. A crucial 
aspect of zero-point fluctuations of the quantum ground state, 
is that they cannot supply energy, but can only contribute to 
processes where energy is absorbed by the mechanics. This is 
completely different from classical noise. Techniques that at- 
tempt to measure zero-point motion without being sensitive to 
this aspect (i.e., standard continuous linear position detection) 
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FIG. 1: Multi-mode Optomechanical resonator, a, A scanning 
electron micrograph of the sihcon nanobeam optomechanical cavity, 
b, and c, show finite-element-method (FEM) numerical simulations 
of the first- and second-order optical modes of the cavity which are 
used for cooling and probing the mechanical motion, respectively, 
d, FEM numerical simulation of the coupled breathing mechanical 
mode. 



can always be interpreted classically and described by some 
effective temperature. 

A more direct method of thermometry and characteriza- 
tion of the quantized nature of a mechanical oscillator, one 
particularly suited to small {n) and utilized in the above- 
mentioned trapped atom experiments (SHSl, is referred to as 
motional sideband spectroscopy. This method relies on the 
fundamental asymmetry in the quantum processes of absorp- 
tion from (proportional to {n)) and emission into (proportional 
to (^) + 1) the mechanical oscillator of phonons. In the case 
of atomic systems, this asymmetry can be measured in the 
motionally-generated Stokes and anti- Stokes sidebands in ei- 
ther the fluorescence or absorption spectrum of the atom. The 
ratio of the Stokes to anti- Stokes motional sideband ampli- 
tudes (((^) + deviates significantly from unity as the 
quantum ground state is reached {{n) 0), and provides a self- 
calibrated reference for the phonon occupancy. In the present 
experiment we cool a nanomechanical resonator to near its 
quantum ground state, and measure the asymmetry in the mo- 
tional sidebands utilizing a form of resolved sideband spec- 
troscopy based upon the filtering properties of a high-Q optical 
cavity with linewidth narrower than the mechanical frequency. 

The cavity optomechanical system studied in this work con- 
sists of a patterned silicon nanobeam which forms an optome- 
chanical crystal (OMC) |30| capable of localizing both opti- 
cal and acoustic waves (see Fig. [T]). The cavity is designed 
to have two optical resonances, one for cooling and one for 
read-out of mechanical motion. The cooling mode is cho- 
sen as the fundamental mode of the patterned nanobeam cav- 
ity, with a frequency 0)^/2:: = 205.3 THz and a correspond- 
ing free- space wavelength of Xc = 1460 nm. The read-out 
mode is the second-order mode of the cavity with co^/27r = 
194.1 THz (kr = 1545 nm). An in-plane mechanical breath- 
ing mode at co^/27r = 3.99 GHz, confined at the center of the 
nanobeam due to acoustic Bragg reflection, couples to both 
optical resonances with zero-point optomechanical coupling 
rates gc/2n = 960 kHz and gr/2n = 430 kHz for the cooling 
and read-out modes, respectively (gc and gr are measured for 



each cavity mode from optically-induced damping of the me- 
chanical mode as in Ref. |i5J). The Hamiltonian of the coupled 
system is given by 

H = h{(Or^grX/Xzpf)a^d^h{(Oc^gcX/Xzpf)c^C 

-^hcOmb^b, (1) 

where we have defined c (c'^) and d (d'^) to be the annihilation 
(creation) operators for photons in the cooling and read-out 
modes of the optical cavity. The optical modes are coupled to 
the mechanical displacement operator of the breathing mode, 
X = Xzpf{b'^ -\-b), where P and b are the creation and annihi- 
lation operators for phonons in the mechanical resonator, and 
Xzpf is the zero-point fluctuation amplitude of the mechanical 
motion. 

An illustration of the experimental apparatus used to cool 
and measure the nanomechanical oscillator is shown in Fig.|2] 
In order to pre-cool the oscillator, the silicon sample is 
mounted inside a Helium flow cryostat. For a cold finger tem- 
perature of 6.3 K, the bath temperature of the mechanical mode 
is measured to be 18 K (corresponding to a thermal phonon 
occupation of Uh = 94 phonons) through calibrated optome- 
chanical thermometry as described in the Methods section and 
Ref.|5|. At this temperature the mechanical damping rate to 
the thermal bath is found to be y//27r = 43 kHz, corresponding 
to a intrinsic quality factor 2^ = 9.2 x 10^. The optical reso- 
nances of the OMC cavity are measured to have total damping 
rates of Kc/27i = 390 MHz and K^/27r = 1.0 GHz for the cool- 
ing and read-out modes, respectively. An optical fiber taper, 
formed from standard single mode optical fiber, is used to op- 
tically probe the OMC cavity via evanescent coupling. The 
component of extrinsic damping resulting from coupling to 
the optical fiber taper waveguide is measured to be Ke^c/^TZ = 
46 MHz for the cooling mode and Kg,r/27r = 300 MHz for the 
read-out mode. 

As alluded to above, resolved sideband cooling in optome- 
chanical cavities follows physics which is formally similar to 
the Raman processes used to cool ions to their motional ground 
state (61 . A cooling laser, with frequency CO/ = 0)^ — 0)^, is 
tuned a mechanical frequency red of the cooling cavity reso- 
nance of the OMC, giving rise to an intra-cavity photon pop- 
ulation He at frequency CO/. The oscillations of the mechanical 
system cause scattering of the intra-cavity cooling beam laser 
light into Stokes and anti-Stokes sidebands at co^ — 2co^ and 
cOc, respectively. Since the anti-Stokes sideband is resonant 
with the cavity at co^, and < co^, the anti-Stokes optical up- 
conversion process is greatly enhanced relative to the Stokes 
down-conversion process. This leads to cooling of the me- 
chanical mode and can be modeled effectively as an additional 
mechanical damping term of Jc = ^g^c^c/^c due to coupling to 
a near-zero temperature bath represented by the cooling laser. 
Assuming a deeply resolved sideband system (k^/co^ <C 1), 
the back-action cooled mechanical mode occupancy is approx- 
imately given by {n)c = ym/iyi^yc) CSIISl. 

Optical scattering of the intra-cavity light field can also be 
used to read out the motion of the coupled mechanical oscil- 
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FIG. 2: Experimental setup. The silicon OMC cavity is mounted into a continuous flow Helium cryostat and optically coupled via an optical 
fiber taper probe. Two narrowband lasers (New Focus Velocity series; linewidth ~ 300 kHz) are used to independently cool and read out the 
motion of the breathing mechanical mode of the OMC cavity. The 1500 nm (read-out) and 1400 nm (cooling) laser beams are passed through 
variable optical attenuators (VOAs) to set the laser power, and combined at an optical fiber wavelength multiplexer (?i— MUX) before being 
sent into the cryostat through an optical fiber. Transmission of the 1500 nm read-out beam through the OMC cavity, collected at the output end 
of the optical fiber, is filtered from the 1400 nm cooling beam light via an optical fiber bandpass filter, pre-amplified by an Erbium-doped fiber 
amplifier (EDFA), and detected on a high-speed photodetector (PD2) connected to a real-time spectrum analyzer (RSA). Other components 
are labeled as amplitude-modulation (a-m) and phase-modulation ((|)-m) electro-optic modulators, fiber polarization controller (FPC), swept 
frequency radio-frequency signal generator (rf-sg), lock-in amplifier (lock-in), optical switches (SW), and calibrated optical power meter 
(power meter). Further details of the experimental apparatus and measurements are in the Methods section. 



later. For a quantum harmonic oscillator, the autocorrelation 
Gxx{t) = {x{t)x{0)) of the position |3 | is given by, 

Gx.(0 Azpf = + {{n) + l)^-'""-'-^!'! , (2) 

(y is the total mechanical damping rate) with its Fourier trans- 
form, the displacement power spectral density (PSD), equal to 

S i^Mx^ - ^(^) 



+ 



y(W + i) 



(co^-co)2 + (y/2)2' 



(3) 



The complex character of Gxx{t) and the resulting asymmetry 
in frequency of S^xiS^) (illustrated in Fig. |3^) arise from the 
non-commutativity of position (x) and momentum {p) opera- 
tors in quantum mechanics and the resulting zero-point fluc- 
tuations. This absorption-emission asymmetry has absolutely 
no classical analogue. Of course, at high phonon occupation 
numbers where {n) ^ (^) + 1, the classically symmetric spec- 
tral density is recovered. Since the optical cavity frequency 
is linearly coupled to the position of the mechanical oscilla- 
tor, the displacement noise spectrum is imprinted on the pho- 
tons leaving the cavity and can be measured optically, with 
the Stokes (anti-Stokes) motional sideband proportional to the 
positive (negative) frequencies of 5'xx(co). 

Consider a read-out laser with frequency CO/^ and detuning 
A = co^ — CO/^ from the read-out cavity mode. The optical power 
spectrum about CO/^ of the transmitted read-out beam leaving 



the cavity is given by fTSl, 



^(co) = |^,outp8(co)- 



2nKr (co^-co)2 + (y/2)2 



27rK, (co^ + co)2 + (y/2)2' 



(4) 



Here A^^ and A^^ are the detuning-dependent anti-Stokes and 
Stokes motional scattering rates, respectively, of the read-out 
laser, given by 



(A±co^)2 + (k,/2)2 



(5) 



As illustrated in Fig. [3}) and c, the optical read-out cavity can 
be used to selectively filter the positive or negative frequency 
components of 5'(co). For a detuning A = — co^ for the read-out 

(r) (r) 

laser. Ay :^A_\ resulting in a Lorentzian signal with area /_ 
proportional to (n) centered at the mechanical frequency 
in the PSD of the photocurrent generated by the transmitted 
read-out laser. Conversely, a detuning of A = co^ results in 

(r) (r) 

A_ AY , producing a signal of area /+ proportional to (n). 
Comparison of the area under the Lorentzian part of the mea- 
sured photocurrent PSD of the transmitted read-out laser for 
detunings A = ±0)^, can then be used to directly infer the me- 
chanical mode occupancy. 
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FIG. 3: Optical filtering of motional sidebands, a, Displacement 
noise PSD, Sxx, of a quantum simple harmonic oscillator, plotted 
against -co for clarity, b. Scheme for measurement of the down- 
converted (Stokes) motional sideband. Here the read-out laser (green 
vertical arrow; frequency CO/^) is detuned a mechanical frequency 
above that of the read-out cavity resonance (green solid curve), c. 
Corresponding scheme for measurement of the up-converted (anti- 
Stokes) motional sideband. The linewidth of the readout cavity (k^) 
and the mechanical resonance (y) are indicated. Insets to b and c show 
a zoomed-out spectra indicating the relative frequency of the cooling 
cavity mode and cooling laser. 



This simple argument neglects the back-action of the read- 
out beam on the mechanical oscillator. In particular, the me- 
chanical damping rate becomes detuning dependent, with y± = 

(7/ +yc)(l ±Cr) for A = ±0)^. Here C, = \A^l^ -A^l^l/iyi^Jc) 
is the effective cooperativity of the read-out beam in the pres- 
ence of the strong cooling beam, and is found from the mea- 
sured spectra by the relation 

C, = ^+~^~ . (7) 

The back-action of the read-out beam also results in a corre- 
sponding change in the phonon occupancy with read-out beam 
detuning, given by {n)± = {n)c/{l ±Cr) for A = ±0)^, where 
{n)c is the mechanical mode occupancy in the presence of only 
the cooling beam. In our measurements this dependence is 
made small by ensuring that the read-out beam is much weaker 
than the cooling beam, so that <C 1 . Adding in a correction 
for the read-out laser back-action, one finds the following rela- 
tion between the measured motional sidebands and the phonon 



occupancy of the cooled mechanical oscillator, 

where for <C 1 we recover the standard relation given in 
Equation d6]). 

Figure [4] summarizes the measurement results of the cali- 
brated mechanical mode thermometry and motional sideband 
asymmetry for the silicon nanobeam OMC cavity. These mea- 
surements were performed with the cooling laser locked a me- 
chanical frequency to the red of the fundamental mode of the 
OMC cavity, with the cooling laser power swept from Uc ^ I 
to 800. For each cooling laser power, the read-out laser is 
used to both estimate the mechanical mode phonon occupancy 
and to compare the motional sideband amplitudes. With the 
read-out laser set to a detuning from the read-out cavity mode 
of A = co^, a Lorentzian spectrum with linewidth y+ and in- 
tegrated area /+ is measured in the read-out laser photocur- 
rent PSD, from which a mode occupancy of (n)^ is inferred 
from independently calibrated system parameters (see Meth- 
ods and Ref.|5|). Similarly, by placing the the read-out laser 
at A = — co^ we obtain spectra with linewidth y_ and inte- 
grated area /_ cx (^)_ + 1, from which we estimate (n)-. FIT 
spectroscopy ll22ll of the read-out optical cavity mode is also 
performed, providing a second, and more accurate (given the 
better signal-to-noise of such measurements) estimate of the 
mechanical damping rates y± . 

From the measured mechanical damping rates, the read-out 
cooperativity Cr is found using relation (|7]), and plotted in 
Figure |4^. The ratios ^ = y_/y+ and {n)^/{n)- are plotted 
in Figure |4]), and found to be approximately equal in accor- 
dance with theory. From {n)±, the estimated phonon occupa- 
tion number in the absence of a read-out beam, {n)c, is calcu- 
lated and plotted in Figure [4]: versus the mechanical damping 
rate y = (y+ +y_)/2. As expected, {n)c drops approximately 
linearly with y, reaching a minimum value of approximately 
2.6 ± 0.2 phonons, limited in this case by available cooling 
laser power. As mentioned, the read-out beam was kept low 
in power so as to minimize the back-action in comparison to 
that of the cooling beam. In particular, for the largest cooling 
powers Cr^2 — 3% and n± are within 10% of {n)c. 

In Figure |4|l, the measured values of the expression r\' are 
plotted versus the calibrated value of {n)c. Also plotted are the 
classical and quantum values of this expression, and l/{n)c, 
respectively. A clear divergence from the classical result of 
r|' = is apparent, agreeing with the deviation due to quan- 
tum zero-point fluctuations of the mechanical oscillator. At the 
largest powers, we measure asymmetries in the motional side- 
band amplitudes on the order of 40% in agreement with the in- 
ferred {n)c = 2.6 phonons from calibrated thermometry. This 
deviation is directly apparent in the measured spectra, shown 
for {n)c = 85, 6.3, and 3.2 phonons in Figure [4^-g, with the 
shaded region corresponding to the noise power contribution 
due to quantum zero-point flucuations. 

While the quantum nature of a mechanical resonator will 
come as little surprise to most physicists, its observation 
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FIG. 4: Measurement of motional sideband asymmetry, a Plot of the cooperativity of the read-out beam as a function of damped mechanical 
linewidth. b Plot of the measured ratios y-/y+ (blue o) and / («)_ (pink o). c Plot of the laser cooled breathing mode phonon occupancy, 
{n)c, as a function of the optically damped mechanical linewidth, y. The dashed line is the predicted phonon number y/fz^/y from an ideal back- 
action cooling model in which the thermal bath occupancy (nfy) and the intrinsic mechanical damping rate (y/) are constant and independent 
of the cooling beam power. In practice both ni, and y/ can vary with cooling beam power due to optical absorption in the silicon nanobeam 
cavity |5 |, which accounts for some of the deviation of the measured phonon occupancy from the ideal cooling curve. Vertical error bars in 
b and c indicate uncertainty in the calibrated phonon occupancy due to uncertainty in the system parameters and a 95% confidence interval 
on the area from the Lorentzian fits, d Plot of the asymmetry (if]') in the measured Stokes and anti-Stokes sidebands of the read-out laser 
for each calibrated measurement of {n)c. The horizontal error bars arise from a 2% uncertainty in the transmitted read-out laser beam power 
between detunings A = ±00^, and a 95% confidence interval in the Lorentzian fitting parameters of the measured spectra from which I± are 
calculated. The vertical error bars corresponding to uncertainty in {n)c are the same as those plotted in c. The classical (blue curve) and 
quantum mechanical (pink curve) relations for the sideband asymmetry, given by and \/ {n)c, respectively, are also plotted, e-g. Plot of the 
measured Stokes (red curve) and anti-Stokes (blue curve) read-out beam spectra for (from top to bottom) {n)c = 85, 6.3, and 3.2 phonons. 
For clarity, we have divided out the read-out back-action from each spectra by multiplying the measured spectra at detunings A = ±C0m by 
y± such that the ratio of the areas of the Stokes and anti-Stokes spectra are given by /_y_//+y+, which approaches unity in the classical 
case. Additionally, we have plotted the horizontal axis in units of y, and rescaled the vertical axis for different {n)c to keep the areas directly 
comparable in each plot. The difference in the Stokes and anti-Stokes spectra, which arises due to the quantum zero-point fluctuation of the 
mechanical system, is shown as a shaded region in each of the plots. 



through the zero-point motion is a significant step towards 
observing and controlling the quantum dynamics of meso- 
scopic mechanical systems. By demonstrating the fundamen- 
tally quantum behaviour of an engineered mechanical nanos- 
tructure, we have shown that realizable optomechanical sys- 
tems have the sensitivity and environmental isolation required 
for such quantum mechanical investigations. 
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METHODS 

Calibrated Thermometry. In this work, an estimate of the mean 
phonon occupation number (n) is accomplished using the read-out 
optical cavity following a calibration procedure outlined in Ref. (5). 
To briefly summarize (see experimental set-up in Fig.[2|, a 1500 nm 
band read-out laser is sent through an amplitude electro-optic mod- 
ulator (EOM), and combined with the cooling beam at a fiber optic 
multiplexer (?i— MUX). PRCs are used to ensure correct polarization 
before the EOMs and the OMC device. Calibration of the detun- 
ing of the read-out (cooling) laser from the read-out (cooling) cavity 
is performed using a form of EIT- spectroscopy f5"22l in which the 
reflection of optical sidebands generated on the read-out laser via a 
swept frequency radio-frequency signal generator (rf-sg) are moni- 
tored on PDl. The intra-cavity photon population generated by either 
the cooling (nc) or read-out («r) laser beam is inferred from measure- 
ments of the power at the input and output of the taper, combined 
with a measurement of the asymmetry in the losses of the taper be- 
fore and after the OMC device. A wavemeter (?i— meter) is used to 
calibrate the optical linewidths of each cavity mode and to lock the 
lasers to within a few MHz of a given frequency. The optomechan- 
ical coupling rate of the read-out cavity, gr, is determined from the 
measured mechanical linewidth versus nr as measured in the trans- 
mitted read-out laser photocurrent spectrum on the RSA. Calibration 
of the mechanical mode occupancy from the read-out laser photocur- 
rent spectrum also requires knowledge of the optical power to detector 
voltage conversion ratio of PD2 (G^) and the amplification furnished 
by the EDFA (Gedfa)- Ge is determined by using SW3 to measure 
the total read-out laser intensity on the calibrated optical power me- 
ter. Gedfa is determined by using SW2 to measure the amplitude of 
a signal generated by the rf-sg, with and without the EDFA in line. 
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